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■ In this talk, I review new developments of QCD calculations of photon and dilepton production 

rates in a Quark-Gluon plasma. All the rates are now known up to 0{as) both for photons and 
dileptons, thanks to the resummation of multiple scatterings. For dileptons, a direct numerical 
^ ' calculation on the lattice attempted recently will also be discussed. 
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1. INTRODUCTION 



> , 

^SJ , Electromagnetic probes (photons or lepton pairs) have production rates that are very sensitive 

■ to the temperature of the medium in which they are produced (the larger the temperature, the 
larger the rate) and are therefore mostly produced during the very early stages of a heavy 



, ion collision. In addition to being produced early, they are weakly coupled to nuclear matter 
I and have therefore a mean free path which is large compared to the typical size of the system 
produced in a collision. This enables them to escape from the system without any reinteraction. 



These two properties combined make electromagnetic signals very good probes of the state of 
the system very early after the initial impact. 

Roughly speaking, a heavy ion collision can be divided into several distinct stages, as illus- 



. trated on figure |T|. Some photons are produced in the initial partonic collisions. These prompt 
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Figure 1. The various stages of a heavy ion collision. 



photons can be calculated using zero temperature perturbative QCD, and they populate the high 
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energy part of the spectrum. Then comes a pre-equihbrium phase, which is usually thought to 
be very short and very poor in quarks and anti-quarks, so that photon production in this phase 
is neglected. In scenarios with a quark-gluon plasma, this is followed by the plasma phase, which 
produces photons at a rate calculable using equilibrium thermal field theory (TFT). Indeed, a 
local equilibrium is usually assumed for the plasma phase, and the local rate calculated using 
TFT is then folded in an hydro dynamical evolution code [ ||, ^ ^ in order to perform the inte- 
gration over space-time. After the confinement phase transition, the system becomes a hadronic 
gas (there can be a mixed phase if the phase transition is first order), for which TFT (now with 
hadronic degrees of freedom) can also be used in order to compute the photon rate. Finally the 
system freezes out and the only photons that are produced afterwards are decay products of 
hadrons. The photon yield observed in detectors is the sum of all these contributions. In the 
rest of this talk, I focus on photons produced in the QGP phase, and on the TFT techniques 
used in order to compute their rate. 

In order to compute the photon production rate in a QGP, one can proceed as follows: make 
the list of the processes contributing up to a given order, add up the corresponding amplitudes, 
square this sum and integrate out all the particles but the photon, properly weighted by the 
appropriate statistical distribution. However, this method becomes rapidly cumbersome since 
one must track by hand the statistical factors, and the interference terms. TFT provides a 
convenient alternative to this approach, which has the advantage of taking care automatically 
of the statistical factors and interferences. In this approach, the photon production rate is 
expressed in terms of the imaginary part of the retarded polarization tensor [ |5|, |6[ : 

"dt^ = (2^)3e^/^'-l • (1) 

The calculation of this retarded imaginary part can be simplified by the use of cutting rules [ 
0, and it is the sum over all the possible cuts that takes care of the interference terms. Note 
that a similar formula exists for lepton pairs, for which one needs to compute the polarization 
tensor of a massive photon (i.e. = — > 0). 



2. HISTORY 

The calculation of thermal photon and dilepton rates has a long and tortuous history. The 
dilepton rate due to the Drell-Yan process (see the diagram on the left of figure ^) was evaluated 
in a QGP in [ ^. Corrections of order 0{as) were considered shortly afterwards, and two 




Figure 2. Real processes contributing to photon and dilepton production up to 0{as). Virtual 
corrections to the first process also contribute at this order. 



problems became apparent. If one tries to calculate the dilepton rate (Q^ > 0) in a plasma 
of massless quarks and gluons, each individual cut contributing to eq. (|l]) contains a mass 
singularity, and it is only after a careful summation of all the real and virtual corrections that 
one gets a finite result [ S, 110, 11 1 . This is nothing but a manifestation of a general property noted 



3 



in 



12|. If one then tries to take the hmit of real photons {Q^ — > O"*") in the above formula, a new 
singularity appears since there are some terms that behave like Imllj-etiuj , q) oc aag ln(u;T/Q^) 



at small Q"^. 

The latter problem was resolved by the resummation of Hard Thermal Loops (HTL). HTL 
are one-loop leading thermal corrections that have the same order of magnitude as their bare 
counterpart when their external momenta are soft (i.e. of 0{^/ctlT)). They were known for quite 
some time in the case of 2-point functions [ |l^ , but a systematic gauge invariant resummation 
was only proposed in 1990 by [ 15, 16|. In the present case, the logarithmic singularity at — > 
is due to the exchange of a soft massless quark. This singularity is screened by the resummation 
of the HTL correction to the quark propagator, which gives the quark a thermal mass m-q of 



order mq ~ 



a.T. 



Taking into account this thermal correction to the quark propagator solves 
the problem and leads to a finite photon polarization tensor [ 17, 18|. For hard photons, it reads: 



Imnret%(a;,q) =4^— r^T 

y 



In 



1 7 (^'('2') 



(2) 



Note that throughout this talk the mass mq is defined to be the asymptotic quark thermal mass, 
irasCfT'^ with Cf = {N^ — l)/2Nc. The numerical factor 5/9 is the sum of the quark 



I.e. m„ 



electric charges squared for 2 flavors (u and d); for 3 flavors (u, d and s), this factor should be 
replaced by 6/9. Regarding the infrared problem, one can see that is replaced by m-q in the 
logarithm as soon as becomes small compared to mq. 

This was thought to be the final answer for the photon and dilepton rates at 0(as), until it 
became clear that some formally higher order processes are in fact strongly enhanced by collinear 



singularities. This was first realized for soft photon production by quark bremsstrahlung [|l^, 2C] 
(left diagram of figure ^). The diagram on the right of figure ^ shares the same property, but 
contributes significantly only to hard photon production [ due to phase-space suppression 
in the case of soft photons. Note that a naive power counting would indicate that these two 
diagrams contribute to O(a^). A common property of these two diagrams is that they have an 
off-shell quark next to the vertex where the photon is emitted, and that the virtuality of this 
quark can become very small if the photon is emitted forward. Again, it is the quark thermal 
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Figure 3. Two processes that are promoted to 0{as) by collinear singularities. 



mass mq that prevents these diagrams from being truly singular. However, contrary to the 0{as) 
diagrams, the singularity is linear instead of logarithmic, and brings a factor T'^/m'^. Combined 
with the that comes from the vertices, these diagrams turn out to be also of order 0{as)- 
In [ and [ |2^, their contribution was derived semi-analytically, with a prefactor that was 

] independently pointed out an erroneous factor 4 in this 



evaluated numerically. [ 22 1 and [ 
numerical prefactor. Finally, it was realized that it can be calculated in closed form 



24, 25 . 
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For 3 colors and 2 light quark flavors, the 0{as) contribution of these two diagrams is exactly: 



ImU.et^ f,{uj,q) = — 

ovr y 



■ rp3 

vr^ — + ujT 

u> 



(3) 



In this formula, the term in 1 /lo dominates for soft photons and comes from the bremsstrahlung 
diagram, while the term in u comes from the second diagram. Because of this term in to, this 
process turns out to dominate the rate of very hard photons {uj ^ T). This was confirmed by 
more realistic evaluations that included this local rate into an hydro dynamical evolution code, 
and there is some speculation that these hard thermal photons could be part of the excess of 
direct photons observed by the WA98 experiment at SPS [ ||, |2|, ^, ^ (see also [ p9|]). 

Note that it is only by accident if this result is so simple. For 3 colors and 3 light flavors, the 
same quantity can still be calculated in closed form (note that the energy dependence is the 
same), but the prefactor is much more involved: 



32 
3^ 



:Li2 



-Lio 



6aas 
9 



vr 



+ ujT 



(4) 



It is also worth mentioning at this point that the purely numerical prefactor is a function of the 
ratio of the quark thermal mass mq to the gluon Debye mass mdebye (the Debye mass quantifies 
the effective range of strong interactions when they are screened by medium effects, and it is 
also of 0{^Joi^T)). In the HTL framework, this ratio is a constant independent of the coupling 
and temperature, that depends only on the number of colors and flavors; for 3 colors and Nf 
flavors, this ratio is mq/mdebye = y^2/(6 + Nf). 



3. LPM EFFECT 

Given the enhancement in the diagrams of figure ^, one may wonder if there are formally 
higher order diagrams that also end up contributing to the same order 0{as)- This was partly 
answered in [ ^] where the resummation of a collisional width F ~ asT\\i{\/ Ug) on the quarks 
in the calculation of the diagrams of figure ^ showed some sensitivity to this parameter at leading 
order, thereby indicating that an infinite series of diagrams must be resummed in order to fully 
determine the 0{as) photon rate. 

In order to explain the issue in more physical terms, it is convenient to define the concept of 
photon formation time. Let me assume that a virtual quark of momentum R = P + Q splits into 
an on-shell quark of momentum P and a photon of momentum Q. The photon formation time 
can be identified with the lifetime of the virtual quark, which is itself related to its virtuality by 
the uncertainty principle. For a small Q^, a simple calculation gives: 



6E = ro 



+ mq 



LO 



(5) 



where the 3-momentum of the photon defines the longitudinal axis. Note that the collinear 
enhancement in the diagrams of figure |3|, due to the small virtuality of the quark that emits the 
photon, can be rephrased by saying that it is due to a large photon formation time. Similarly, 
the sensitivity to the collisional width found in [ ^] occurs if the photon formation time is of 
the same order or larger than the quark mean free path between two soft collisions (this mean 
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free path A is the inverse of the width T). This phenomenon is nothing but a manifestation of 
the Landau Pomeranchuk Migdal (LPM) effect [ ||, ||, |3|. 



The precise nature of the multiple scattering diagrams that must be resummed depends in fact 
on the range of the interactions in the medium. Indeed, if the range i of the interactions is much 
shorter than the mean free path, it is easy to check that only ladder topologies are important, 
in which all the successive scatterings are independent of one another, as illustrated in figure 
^. Indeed, the condition £ <^ X suppresses diagrams with crossed gluons. On the contrary, if 



Figure 4. A ladder correction to bremsstrahlung. 



there are long range interactions in the system, for which ^ ^ A, then arbitrarily complicated 
topologies can a priori contribute. It was found in [ 34| that if one considers contributions to the 



photon polarization tensor topology by topology, then there can be a sensitivity to interaction 
ranges as long as the magnetic scale l/a^T (for a review of the relevant scales and associated 
physics in a QGP, see [ ^), which would render the problem practically intractable. 

A considerable progress was made recently in [ |3^, in which it was shown that there are 
infrared cancellations between diagrams of different topologies, and that these cancellations 
remove any sensitivity to the magnetic scale. Physically, this cancellation can be interpreted as 
the fact that ultrasoft scatterings are not efficient in order to induce the production of a photon. 
As a consequence, only the ladder family of diagrams needs to be resummed in order to obtain 
the complete leading 0{as) photon rate. The resummation of this series of diagram can then 
be performed in two steps summarized in figure 0. The first one is a Dyson equation for the 



Figure 5. Resummation of ladder diagrams. 



photon polarization tensor, whose explicit form is 



imn 



ret 



(Q) ~ aNc 



+ 00 



dpo [np{ro) 



npiPo)] 



3|,|33, H: 



Re 



(6) 



(poro)2 J (27r)2 

with ro = po + 90) "iT-piPo) = ^/i^w{Po/T) + 1) the Fermi-Dirac statistical weight, and where 
the dimensionless function f{p±) denotes the resummed vertex between the quark line and the 
transverse modes of the photon (this is represented by the shaded vertex in the above pictures). 
In the Dyson equation, this function is dotted into a bare vertex, which is proportional to 
The second equation, that determines the value of f{p±), is a Bethe-Salpeter equation that 
resums all the ladder corrections [ 



37 



dH 



(2vr) 



c(/±)[/(p^ + U-/(p±)] 



(7) 
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where tp is the time defined in eq. (0) and where the colhsion kernel has the following expression: 
C{1±) = Tn'^g^yg/l±{l± + "ijgjjyg) [ gj]. Note that in the Dyson equation, the quark propagators 
should be dressed in a way compatible with the resummation performed for the vertex, in order 
to preserve the gauge invariance. It is this dressing on the quark propagators which is responsible 
for the term —f{pj_) under the integral in eq. (^). From this integral equation, it is easy to 
see that each extra rung in the ladder contributes a correction of order asTporQ/Lorn^, in which 
the as drops out. Therefore, all these corrections contribute to 0{as) to the photon rate. Note 
again that the only parameters of the QGP that enter in this equation are the quark thermal 
mass niq and the Debye screening mass rridebye- This integral equation was solved numerically 
in [ and the results are displayed in figure ^. In this plot, 'LPM' denotes the contribution of 
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Figure 6. 0{as) contributions to the photon production rate in a QGP. The parameters used 
in this plot are as = 0.3, 3 colors, 2 flavors and T = 1 GeV. 



all the multiple scattering diagrams, while '2^2' denotes the processes of figure |2|. The single 
scattering diagrams (figure ^) are also given so that one can appreciate the suppression due to 
the LPM effect. 

4. DILEPTON PRODUCTION 

Dilepton production basically suffers from the same problems, and the solution follows the 
same path. Two differences are worth mentioning here. First of all, the Drell-Yan process qq — > 
7* I'^l" contributes if > Arn^. In addition, virtual photons have a physical longitudinal 
mode that contributes to the rate of lepton pairs. The Drell-Yan process has been evaluated in 



[ g] , and the 2 — > 2 processes have been evaluated in [ 

For photon invariant masses of order Q ~ OsT"^ or smaller, one expects also important 
contributions from multiple scattering diagrams. One must now keep track of the non-zero 
Q^, and include also the contribution of the photon longitudinal mode. This is easily done by 
performing the following substitution in eq. (^) [ |39| : 

f^Px • /(PJ - f^P± ■ fiP±) + -7^%9iP^) , (8) 
(Po^o) [poror V\Po'^o\ % 



where the function g{pj_), which describes the coupling between the quark line and the longitu- 
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dinal photon, obeys an integral equation similar to eq. [39]: 



—gip^) = 2^\poro\ + ATrosCfT 



(9) 



Note that the contribution of the longitudinal mode of the photon vanishes trivially when 
0, as it should. This new integral equation can also be solved numerically, and the resulting 
dilepton rate (for the same parameters as in figure ^ and a total energy of the pair set to 
(jj = 5 GeV) is plotted in figure ^. One can see that the multiple scattering corrections are 

a,=0.3, 3 colors, 2 flavors, T=lGeV, co=5GeV 
Threshold = 2.3 GeV 
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Figure 7. 0{as) contributions to the dilepton production rate in a QGP. 



important for all pair masses below the threshold of the Drell-Yan process. Note also that the 
threshold of the tree-level process is completely washed out when multiple rescatterings are 
resummed. 

5. QUASIPARTICLE MODELS 

We have emphasized several times the fact that the only properties of the QGP that these 
rates depend on are the quark thermal mass and the Debye mass ?n.debye- So far, these two 
masses have been taken in the HTL approximation, for which the ratio of the two masses is 
independent of T and as- However, simple arguments indicate that this ratio cannot remain 
constant when the mass ruq becomes large, which may happen at moderate temperatures for 
which the coupling constant is rather large. Indeed, the Debye screening is due to the possibility 
for a test charge to polarize the medium surrounding it in order to screen its charge. This 
process becomes difficult to achieve when the quasiparticles in the medium become very heavy, 
and for this reason mdebye should become very small if niq increases. This is indeed what one 
finds by calculating the Debye mass at 1-loop, with massive particles running in the loop. 

In practice, one could obtain the mass niq from a quasiparticle fit to the lattice entropy, as 
has been done in [ This is illustrated in figure ^. On the left plot, one can see that this 
fit reproduces perfectly the entropy measured on the lattice, even for temperatures very close 
to the critical temperature Tc- The quasiparticle mass (here it is a gluon mass nig, but this is 
proportional to niq) needed for that fit is shown by a dashed line in the right plot, and one can 
see that it becomes very large near Tc- The solid line is the Debye mass calculated at 1-loop, 
with the mass obtained from the previous fit used for the particle in the loop. The dots are 



8 



0.75 - 




0.25 



Figure 8. Extraction of niq from a fit of the lattice entropy, and comparison of the resulting 
TTT-debye with the Debye mass measured on the lattice (see the text for more details). 



the values of the Debye mass measured on the lattice, and one can see that at least the trend 
is remarkably well predicted by this very simple model. Since photon rates often need to be 
evaluated at temperatures that are not very large compared to Tc, it could be important to take 
the values of mq and rridebye from this model rather than from the HTL approximation. 

6. LATTICE CALCULATIONS 

Recently appeared the first attempt to calculate directly on the lattice the production rate 
of dileptons in a quark- gluon plasma. In fact, the principle of this calculation has been known 
for a long time: one should start from the Euclidean correlator of two vector currents n(T, x) = 
{jfj,{0,0)j'^{T,x)), where r G [0,1/r] is the Euclidean time. Next, one obtains n(T, q) by a 
Fourier transformation of the spatial coordinates, and the imaginary part of the real time self- 
energy is then related to this object by a simple spectral representation: 

X /"~, X ^ „ / , cosh(w(r - l/2r)) 
Jo smh(cj/2r) 

In fact, this equation uniquely defines lm.Il-ret^ ^{u , q) if n(r, q) is known for all r G [0, 1/T] and 
if one prescribes the behavior of the solution at large uj. 

The main problem in lattice calculations is that one knows the function n(r, q) only on the 
discrete points of the lattice, which prevents from determining uniquely the solution. This 
problem has been reconsidered recently using the Maximum Entropy Method [ which 
is a way to take into account prior knowledge about the solution (positivity, behavior at the 
origin, etc..) in order to determine the most probable solution compatible with the lattice data 
and with this a priori information. The result obtained for static dileptons {q = 0) via this 
method is displayed in the figure ^, for two different values of the temperature. Note that this 
is a quenched lattice simulation. This result displays several interesting properties. At energies 
above AT, the full rate is very close to the contribution of the Born term, while at energies 
smaller than 2>T it drops to extremely small values. In addition, when plotted against w/T, the 
curves for the two temperatures fall almost on top of one another, indicating that the result 
scales like a universal function of cj/T, at least within the errors. 

In fact, the suppression at small u: has attracted a lot of interest because it is not what 
one would expect from perturbation theory: the resummation of thermal masses would indeed 
produce a drop of the Born term because of threshold effects, but there are some higher order 
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Lattice calculation for static lepton pairs 
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Figure 9. Lattice result for the production rate of static dileptons. 



processes that do not have a threshold and that should fill the spectrum at small w. As of 
now, there are arguments indicating that both the perturbative calculations and the lattice 
calculation are incorrect at small uj. If one evaluates eq. (10) at r = l/2r, one gets a sum 
rule: dojluill^^t^ ^{uj , q)/sinh(a;/2T) = n(l/2r, q) < oo, which is violated by all the existing 
perturbative calculations (they give an infinite result) because none of them includes the strong 
dissipative effects that appear when one enters in the hydro dynamical regime (a; ^ 0). 

On the other hand [ 43|, from the electric conductivity: cJei = limtj_>oImni.et(w,0)/6a;, one 
obtains Imnret('-^, 0) oc w when u; — > 0. This implies that the static dilepton rate should diverge 
when a; ^ 0. Unless the electric conductivity in quenched QCD is nearly zero for some (yet to be 
explained) reason, the lattice dilepton rate disagrees with this prediction at small w. Note that 
'small' in these considerations means an iv small enough to be in the hydrodynamical regime, 
i.e. UJ < g^T. In a strong coupling theory, this regime could start as early as w ~ T. 



7. CONCLUSIONS 

The full (D{as) photon and dilepton rates have now been calculated using thermal QCD. This 
required to resum all the diagrams involving multiple rescatterings as it turns out that the LPM 
suppression plays a role at this order. A possible improvement of this perturbative calculation 
could come via the use of quasiparticle models adjusted to reproduce thermodynamical quantities 
determined on the lattice, in order to get more realistic values for the various mass parameters 
that describe the internals of the QGP. 

Another recent development is the direct lattice evaluation of the static dilepton rate. More 
work is still needed in this area in order to fully understand the discrepancy with the perturbative 
approach at low energy, and also to extend this program to real photons. 

Finally, one should also mention that the way the pre-equilibrium phase is treated (or, more 
accurately, ignored) is probably not correct as there are indications that the kinetic equilibration 
time might be as long as a few fermis [|4^. Determining how quarks are dynamically generated 
from the initial gluons is a topic that certainly deserves much more attention, and may influence 
strongly the photon yields on obtains at the end. 
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